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ABSTRACT 


The capabilities provided by the PATRAN-COSMIC/NASTRAN interface are 
discussed. While the translator capabilities give some indication of the 
interface quality between the two proqrams, certainly there are otner 
attributes to be considered. The ideal interface would he a user-transparent 
union of the two programs so that the engineer could move from one program to 
the other fluently ar.d naturally. Hence, a valid assessment of the interface 
completeness must consider how close the current capabilites are to the 
idealized case. An example problem is presented to demonstrate how 
COSMIC/NASTRAN and PATRAN can be used together to meet the requirements of an 
actual engineering application. 


INTRODUCTION 

PATRAN [1] is a three-dimensional solid modeling and finite element pre- 
and postprocessing program developed and marketed by PDA Engineering in Santa 
Ana, California. Using the latest interactive computer graphics technology, 
PATRAN provides a visual means to define a finite element model and its 
environment, and review its resultant model behavior. Although PATRAN offers 
linear statics and eigenvalue analysis capability as an option, a more 
detailed analysis is often required. More than 25 translator programs have 
oeen developed that move data among PATRAN and external finite element, finite 
difference, and boundary element analysis programs [2]. It is these 
interfaces that bring the advanced analysis methods available in 
general-purpose finite element programs to the PATRAN modeler, and one 
important interface is to COSMIC/NASTRAN [3]. 

The COSMIC/NASTRAN interface is comprised of two translators: an 
analysis translator and a results translator. The analysis translator accepts 
a PATRAN neutral output file and produces a NASTRAN bulk data file. The 
forward translator is named PATCOS to indicate the translation from PAT RAN to 
COSM IC/NASTRAN. For postprocessing within the PATRAN system, it is necessary 
to execute a second program after the analysis. This inverse translator is 
called COSPAT, because it accepts the COSMIC/NASTRAN results and produces 
files in a format which PAT RAN can read. CdSPAT also provides the option of 
translating an existing analysis input deck into a PATRAN neutral file. Both 
analysis and results translators are easy to use in that they are menu- driven, 
execute quickly, and are highly automated. The PATRAN-COSMIC/NASTP.AN interface 
is illustrated lr, Figure 1. 
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ANALYSIS MODELING 


PATCOS can produce 59 different COSMIC/NASTRAN bulk data card types, 
including 29 different finite elements. The list of bulk data cards supported 
is presented in Table 1. The complete analysis model can be defined in PATRAN 
to include not only the finite element mesh but also element properties, 
boundary conditions end applied loads. PATRAN also provides visual model 
verification techniques to check parameters such as element aspect ratio or 
property assignment. 

There are several ways to add the required EXECUTIVE and CASE CONTROL 
cards to the bu’k data deck. The simplest method is to use a text editor and 
either (1) insert the control cards manually, or (2) copy them in from an 
existing file. 


RESULTS RECOVERY 


The results (or inverse) translator, COSPAT, provides the capability to 
reformat analysis results into PATRAN compatible files. The input file to 
COSPAT is a binary 0UTPUT2 file generated during COSMIC/NASTRAN execution. The 
output files contain nodal displacements, element centroidal and nodal stress, 
strain or temperature. To generate an 0UTPUT2 file during a COSMIC/NASTRAN 
execution of a static analysis (SOL 1), the following directives must be added 
to the COSMIC/NASTRAN input deck. Similiar ALTER sequences can be used for 
other solution approaches: e.g. 

ALTER 143 $ 

0UTPUT2 0UGV1 ,0ES1//C,N,-1/C,N,11/V,N,Z $ 

ENDALTER 

CEND 

DISP = ALL 
STRESS = ALL 

PATRAN can then be used to combine the results information with the 
finite element model to evaluate model behavior. Using graphical methods such 
as animation and color bring unparalleled insight into understanding the 
voluminous output data often generated from a complex analysis. The engineer 
can interactively view any piece of the model from any perspective and see 
color variations of results. Hence one simple exercise for PATRAN might be to 
crack open an egg and display the variation in circumferential stress on the 
inside of the shel 1 . 
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MODEL TRANSLATION 


COSPAT also provides the capability to reformat COSMIC/NASTRAN analysis 
input data into PATRAN neutral format. This method not only brings all PATRAN 
postprocessing capabilities to bear on existing COSMIC/NASTRAN models 
(generated by other means o* by hand), but also suggests a medium in which to 
convert a COSMIC/NASTRAN input deck to that of another analysis program 
format. The COSPAT model translator currently recognizes the 58 different 
card types as listed in Table 2. 


AN APPLICATION 


While the translator capabilities give some indication of the interface 
quality between the two programs, it is clear that there are other attributes 
to be considered. The ideal interface would be a user-transparent union of the 
two programs so that the engineer could move from one program to the other 
fluently and naturally. Hence, a valid assessment of the interface 
completeness must consider how close the current capabilites are to the 
idealized case. This can be explored through the use of an example problem. 

One common step in the analysis of a part or structure is to determine 
the structural response to an imposed thermal environment. This involves 
creating a thermal model subjected to various thermal conditions, solving for 
a resultant temperature distribution, and then calculating the stresses 
induceo by the applied thermal loads. Both a thermal and a structural analysis 
must be performed, very often using a structural analytical model of a 
different mesh density than that of the thermal analysis. The following 
example demonstrates how PATRAN and COSMIC/NASTRAN can be used together to 
model, analyze, and evaluate a thermal stress problem. 


Problem Definition 

The example chosen is a pipe with stepped cooling fins subjected to an 
internal fluid temperature of 500"F and an outside ambient temperature of 70*F 
(Figure 2). The structure consists of a 6-inch (0.5 ft) stainless steel pipe 
with 14.4-inch (1.2 ft) cooling fins. The emphasis here is not on the actual 
dimensions and properties of the model but rather the technique applied. The 
objective is to use axisymmetric elements with heat boundary elements to solve 
for a linear steady-state temperature distribution. Then, determine the stress 
induced by the resulting temperatures on two different structural models: (1) 
a finite element mesh matching the thermal mesh, and (2) a non-uniform 
structural mesh with higher element density. 
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Thermal Model 


The thermal model is represented with axisymmetr ic trapezoidal ring 
elements (CTRAPRG). The only material property required for this element is a 
thermal conductivity equal to 15 Btu/hr-ft- *F. The convection surfaces are 
defined by heat boundary line elements (CHBDY). The material property for this 
element is the convective film coefficient equal to 12 Btu/hr-ft-*F on the 
inside of the pipe and 3 Btu/hr-ft-^F on the outside of the fin. The 
appropriate PATRAN directives are used to define these data (Table 3), namely 
element connectivity (CFEG), element geometric properties (PFEG), and material 
properties (PMAT). The only other value to be defined is the ambient 
temperatures referenced by the convection surfaces. Because PATRAN originated 
from structural analysis, there is currently no straightforward way to define 
this relationship. Hence, a special technique must be used. This technique 
consists of first constraining scalar points to the appropriate fixed 
temperatures and then 'eferencing these scalar points in the element associate 
data field as the convection surface elements are created. The complete list 
of PATRAN directives required to generate the thermal model is listed in Table 
3. The final model is illustrated in Figure 3. 

PATCOS is then executed to translate the PATRAN neutral output file into 
COSMIC/NASTRAN bulk data. A minor coding modification translated the element 
data associated with the heat boundary elements as the ambient reference 
points. The EXECUTIVE and CASE CONTROL cards were added and included two 
important directives: (1) THERMAL(PUNCH) = ALL to write the resultant nodal 
temperatures to a punch file and (2) 0UTPUT2 to save the temperatures for 
PATRAN postprocessing displays. A sampling of the COSMIC/NASTRAN input file is 
shown in Table 4. 

COSPAT can be executed to read the 0UTPUT2 file and generate PATRAN 
postprocessing files. Figure 4 illustrates the temperature distribution that 
resulted from the linear steady-state heat transfer analysis. 


Structural Model (Mesh 1) 

The first structural model is created by slightly modifying the thermal 
model using PATRAN. The heat boundary elements are deleted, the structural 
properties for the CTRAPRG elements are defined, and axisymmetric boundary 
conditions are added. For this case, Young's modulus is given as 194,400 
lb/ft, Poisson's ratio is 0.3, and the coefficient of thermal expansion is 
0.00009 ft/ft- F. 

PATCOS is then used to generate the corresponding COSMIC/NASTRAN bulk 
data. This time the appropriate EXECUTIVE and CASE CONTROL cards are added, as 
well as a single point constraint card to eliminate all degrees of freedom 
except the in-plane axisymmetric motion. Also, because the structural mesh 
corresponds identically to the preceding thermal mesh, the punch file 
generated during the thermal analysis can be appended directly to the 
structural COSMIC/NASTRAN input file. Clearly, this is a most desirable 
situation and, for this reason, many thermal models are commonly over-modeled 
such that a structural analysis can be performed easily at a later time. 
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The 0UTPUT2 file generated during the COSMIC/NASTRAN run is processed by 
COSPAT and PATRAN is used to plot the deformed shape (Figure 5) and stress 
contours (Figure 6). 


Structural Model (Mesh 2) 

The next step is perhaps the most difficult, in that a new structural 
mesh is to be defined but subjected to the temperature distribution calculated 
in the inital thermal analysis. Many schemes have been devised to interpolate 
resultant temperature fields onto new structural mesnes. The inherent 
capability to the PATRAN approach of defining a model's environment in terms 
of the model geometry and not the finite element mesh shows great potential to 
solving this problem. However, direct methods have not yet been developed. 

The approach taken here is to define PATRAN data patches (surfaces) that 
lie in the thermal distribution field. PATRPN data entities (lines, surfaces, 
and volumes) can be created to represent a parametric cubic variation of any 
scalar function [4]. These data surfaces can then be applied to the geometry 
independent of the finite element mesh density. This is made possible by 
taking advantage of the axisymmetric nature of the problem. It is noted that 
the spatial orientation of the model is completly defined by X- and 
Z-coordinates and that the Y-coordinate is always zero. Hence, a 100-line 
FORTRAN program, developed in two hours, combines the spatial model 
coordinates with the nodal temperatures to define PATRAN geometric points 
(grids) in a PATRAN neutral file. The three coordinates of these grid points 
are specified as X, TEMP, Z so that they actually lie in the temperature 
field. 


The next step is to delete all finite elements and nodes in the PATRAN 
axisymmetric model, leaving only the model geometry. The temperature neutral 
file is then combined th the PATRAN model and these new grid points are used 
to define PA T RAN data patches. The cubic data patches are defined by 
specifying scalar values that lie in tre distribution field. Hence, it is an 
easy task to have PATRAN extract the f-coordinates of selected grid points to 
define the desired data patch (scalar function). 

A new structural mesh is defined with a non-uniform element density 
(Figure 7). Then, each patch is loaded with its corresponding data patch and 
each finite element node is assigned automatically a temperature depending on 
its location in the data patch function (Figure 8). Geometric patches were 
constructed to show the temperature distribution. 

PATCGS then translates this new model and loading combination into 
COSMIC/NASTRAN bulk data. The EXECUTIVE and CASE CONTROL cards are added and 
the analysis is run. COSPAT is used to generate the postprocessing files that 
lead to the final Z-X shear stress plots for the critical area shown in Figure 
9. 
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CONCLUDING REMARKS 


An overview of the current capabilities provided by the 
PATRAN-COSMIC/NASTRAN interface is described. While the extent of these 
capabilites is important in assessing the interface quality, the example 
problem illustrates the wide range of flexibilty provided. PDA Engineering, in 
its short history of interface development, has made major strides in linking 
up the unique analysis tools found in many different programs. Improvements 
have been made which were primarily due to a close working relationship with 
the engineering community. We were thus able to define, and meet, the risinq 
needs and expectations of the analysts. 

The future holds a great promise for more efficient, complete, and 
sophisticated methods in software and hardware networks, : o solve a growing 
diversity of engineering problems. PDA’s interface developments already are 
cone Spring new approaches to thermal modeling, substructuring, and composite 
material modeling. It is clear that PDA is committed to providing the 
engineering community with the highest quality modeling and analysis software 
through more efficient, complete, and accurate data definition and transfer -- 
very often based on user recommendations. 
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Table 1 

COSMIC/NASTRAN Card Types Supported By PATCOS | 59) 


Coordinate Frame' 

Element Properties 

Node Forces 

C0RD2C 

PBAR 

FORCE 

C0RD2R 

PHBDY 

MOMENT 

C0RD2S 

PIHEX 



PQDMEM1 

Specified 


PQDMEM2 

PQDPLT 

PQUAD1 

Node Displacements 

Node Coordinates 

PQUAD2 

SPC 

GRID 

PSHEAR 



PTRBSC 

Constraints 


PRTIA1 


Element Definitions 

PTRIA2 



PTRIM6 

SPC1 

CBAR 

PTRMEM 


CELAS2 

PTRPLT 


CHBDY 

CHEXA1 

CHEXA2 

PTRSHL 

Temperatures 

CIHEX1 


TEMP 

CIHEX2 

CIHEX3 

C0NM2 

Material Properties 


CQDMEM1 

MAT1 


CQDMEM2 

MAT2 


CQDPLT 

MAT3 

Bar Deformation 

CQUAD1 

MAT4 


CQUAD2 

MAT 5 


CROD 

CSHEAR 

CTETRA 

CTRAPRG 

CTRBSC 


DEFORM 

CTRIA1 

CTRIA2 

CTRIARG 

Pressure Loads 


CTRIM6 

PLOAD 


CTRMEM 

PL0AD2 


CTRPLT 

CTRSHL 

CWEDGE 

CNGRNT 

PLOA03 
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Table 2 

COSMIC/NASTRAN Cards Supported By COSPAT (58) 

A. Coordinate frames 
C0R02C, C0RD2R, C0RD2S 

B. Grid points 
GRID .GRDSET 

C. Elements 

BAROR , TRIAi, CBAR . CTRIA3, CROD , CQDMEM1 
CTRBSC, CTRIfRG, CTRIA2, CTRHEM, CQUAOI, CQDMEH2 

CSHEAR, CTRAPRG, CQDMEM, CQUAD2, CHEXA2 , CWEOGE 

CIHEX2, CTRPLT , CTETRA, CQDPLT, CTRIM6, 

CTRSHL, CTRPLT1, CHEXA1, CIHEX1, 

D. Properties 

PBAR , PBEAM , PROD , PTRIA1, PTRMEH 
PQUAD1 , PQUAD2, PQDHEM, PTRIA2, PHEX 
PSOLID, PHBOY , PSHEAR, PSHELL 

E. Materials 

MAT1 , MAT2 , MAT3 , MAT4 , MAT5 , MAT9 

F. Forces and Constraints 

SPC , SPC1 , FORCE , MOMENT, TEMP , DEFORM 
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Table 3 

PATRAN Input Directives To Create Thermal Model 


GO 

1 

1 

A 

2 

VI 

1 

■?0 

GR.1.,.45 
GR,2,,.5f- 
GR,3, , .8*> 

GR,4,,1.75 
LI ,3#,2G, ,1T3,2T4 
LI.#,TR,//.3,1 
LI ,#,TR,// .15,2 
LI ,#,TR,// .05,3 
PA,i#,2L,,lT3,4T6 
GF,P1, ,4/7 
GF.P2, ,4/4 
GF,P3, ,13/2 
CF, PIT#, QUAD/4/7 
LI ,4#,2G, ,6/8/10/5,7/9/4/1 
GF,10L,,7 
GF.4L/7/5, ,4 
GF,8L, ,3 
GF,6L, ,13 
GF,9L,,2 
COLOR, BAR, BLUE 
GRID, 100,, 100 
GRID, 200,, 200 
CF.10L, BAR/2/5 
1 
G 

100 

4 

CF.4LT9, BAR/2/5 
1 
G 

200 

4 

E 

2 

3 


! Initialize PATRAN 

! Change view to X - Z plane 
! Begin geometry creation 


! Finish geometry creation 
! Define thermal mesh 


! Define axisymmetric quad, elements 
! Create geometry for convection boundary 
! Mesh boundary 


! Define scalar points 
! Define convection elements 


! Reference ambient source 


! Use menu sytem to optimize bandwidth 
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Table 3 (Continued) 


N 

2 

1 

Y 

7 

3 

3 

1 

NODE, 100, ADO 
1 
1 

100.0 

5 

NODE, 200, ADD 
1 
1 

200.0 

5 

PISP,N100,AD0,1 
500. C 


! Define s.alar points for ambient temperature 


! Specify ambient temperatures 


DISP,N200,ADD,1 

70.0 

PF.P1T3, QUAD/4/7,, HI ! Assign material types 

PF,10L,BAR/2/5,2,2 

PF,4LT9,BAR/2/5,3,3 

PMAT,1,TIS,15 ! Specify material properties 

PMAT,2,TIS,12 

PMAT,3,TIS,3 

E 

5 ! Create PAT RAN neutral file 

1 
1 

AXISYMMETRIC COOLING FIN (THERMAL MODEL) 

N 

STOP 
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Table 4 

Sampling Of NASTRAN Thermal Analysis Input 

10 DHL, NASTRAN 
APP HEAT 
SOL 1 
TIME 90 
OIAG 14 
ALTER 109 5 

PARAM ' i C,N,N0P/V,N,TRUE=-1 

EQUIV H0UGVI,0UGV1,/TRUE % 

0UTPUT2 ,,,,//C,N,-l/C,H,92/V,N,Z % 

0l)TPU T 2 0UGV1 , ,, ,//C,N,0/C,N,92/V,N,Z $ 

0UTPUT2 ,,,,//C,N,-9/C,N,92/V,N,Z S 

ENDALTER 

CEND 

TITLE = LINEAR STEADY STATE CONDITION THROUGH A COOLING FIN 
SUBTITLE ■ AXISYMMETRIC RING ELEMENTS, FILM HEAT TRANSFER 
OUTPUT 

THERHAL(PRINT, PUNCH) =ALL 
OLOAD = ALL 
SPC = 1 
BEGIN BULK 

JAXISYItlETRIC COOLING FIN (THERMAL MODEL) 

S--BULK DATA CARDS PRODUCEO BY "PATCOS" VERSION 1.6A: 09-JAN-85 15.54:40 


SPC 

1 

100 1 

500 

.000 






SPC 

1 

200 1 

70. 

0000 






GRID 

1 

0.45000 

0 . 


0 . 





GRID 

c 

0.48333 

0 . 


0 . 





GRID 

3 

0.48333 

0 . 


0.05000 





GRID 

4 

0.45000 

0 . 


0.05000 





GRID 

100 

100.0000 

0 . 


0 . 





GRID 

200 

200.0000 

0 . 


0 . 





CTRAPRG 

1 

1 2 


3 

4 

0.000 

1 



CTRAPRG 

2 

2 5 


6 

3 

0.000 

1 



CTRAPRG 

3 

5 7 


8 

6 

0.000 

1 



CTRAPRG 

4 

4 3 


9 

10 

0.000 

1 



CHBDY 

40 

2REV 


26 

22 



E 

40 

+E 40 

100 

100 








CKBDY 

41 

2REV 


22 

18 



E 

41 

+E 41 

100 

100 








CHBDY 

42 

2REV 


18 

14 



E 

42 

+E 42 

100 

100 








CHBDY 

43 

2 REV 


14 

10 



E 

43 

+E 43 

100 

100 








CHBDY 

69 

3REV 


64 

63 



E 

69 

+E 69 

200 

200 








PHBDY 

2 

2 








PHBDY 

3 

3 








MAT4 

1 

15.0000 








MAT4 

2 

12.0000 








MAT4 

3 

3.00000 









ENDDATA 
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Figure 1 

PATRAN - To - COSMIC/ NASTRAN Interface 
(August 1984) 




Figure 2 - Pipe With Cooling Fins 
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Figure 3 - The Thermal Model (also Structural Mesh 1) 











Figure 7 - Structural Model (Mesh 2) 
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